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ABSTRACT 

In this paper we analyse the ahgnment of galaxy groups with the surrounding large 
scale structure traced by spectroscopic galaxies from the Sloan Digital Sky Survey 
Data Release 7. We characterise these alignments by me ans of an extensi on of the 
classical two-point cross-correlation function, developed by iPaz et all ()2008[ ) . We find 
a strong alignment signal between the projected major axis of group shapes and the 
surrounding galaxy distribution up to scales of 30 Mpc h~^. This observed anisotropy 
signal becomes larger as the galaxy group mass increases, in excellent agreement with 
the corresponding predicted alignment obtained from mock catalogues and ACDM 
cosmological simulations. These measurements provide new direct evidence of the 
adequacy of the gravitational instability picture to describe the large-scale structure 
formation of our Universe. 

Key words: galaxies: groups: general, dark matter, large-scale structure of universe. 



1 INTRODUCTION 

The majority of studies of clusters of galaxies have as- 
sumed they are the visible component of a gravitation- 
ally bound system dominated by a dark matter halo. It 
is well known that these objects exhibit mildly aspherical 
shapes, with a slight preference toward prolate forms. More- 
over, it has been shown that their orientations are related 
to the surrou nding structures suc h as filaments and large- 
scale walls llColberg et al.l l2005l: jKasun fc Evrardl 1 2005 : 



Basilakos et al.ll2006l: lAUgood et ani200d: lAltay et al.ll200e 



Aragon-Calvo ct al.' 2007*: 'Bru nino et all BoOTI : iBett et all 
2007; Zhang ct al. 2009). Results from numerica l simula- 
tions by Ivan Haarlem fc van de Weveaertl (| 19931 . see also 
ISplinter et al.lll997h have shown that the origin for such 
alignments comes from the rearrangement of dark matter 
halo axes in the direction from where matter was predom- 
inantly accreted (i.e. last major merger event). Using nu- 
merical simulations, several authors have found that dark 
matter halos tend to be more prola te and aspherical when 
larger halo ma sses are considered (|Kasun fc Evrardl [2OO5I : 
|Paz et al]|2006l . and references therein). This trend is con- 
sistent with the picture described above, and could be eas- 
ily understood as the result of the accretion process driven 
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by mergers. After a merger occurs, dynamical relaxation 
will tend to drive a halo closer to isotropy. Therefore, it 
is expected that dynamically younger systems will be more 
strongly ellipsoidal than dynamically older ones. Given that 
the dynamical age of a system depends on how much time 
has elapsed from the moment of the last halo merger, it is 
expected that high-mass halos will tend to be more elon- 
gated than lower mass halosFl lPaz et al.l (|2006l ) also found 
that this predicted trend is consistent with observational 
results obtained from samples of galaxy groups. More re- 
cently, [Robotha^^eF^al] l|2008h confirms this agreement be- 
tween gal axy group shapes and simulate d dark matter halos. 
Moreover iRagone-Figueroa et al.l (|2010l ) found that less con- 
centrated groups, corresponding to later formation times, 
show lower velocity dispersions and higher elongations than 
groups of the same mass with higher values of concentration, 
which correspond to earlier formation times. Therefore, in 
a statistical sense, the shapes of galaxy groups and their 
alignments could encode information about the formation 
of l arge-scale struct ure (LSS). 

iBinggelil (|1982| ) was the first to investigate alignments 
of pairs of clusters close to each other. Taking 44 Abell 
groups, Binggeli found that galaxies separated by up to 



^ Small halos may start the ir development later than higher mass 
systems l lLagos et al.|[2009l) . but the former will still be s ubject 
to a lower number of major mergers l lDe Lucia et al.ll2006h . 



2 D. J. Paz et al. 



30 Mpc show a strong alignment, while the orientation 
of a given cluster is related to the spatial distribution 
of neighbour syste r ns. W ith some exceptions , for e xample 



Struble & PeeblesI (|l985h or iRhee & Katgert (Il987), most 


papers in the literature confirm this result dde Theiie et al. 


1995 


: IWesd Il989al: iRhee et alJ Il992l: lOnuora & Thomas 


200C 


: IWest et all Il995l). Other authors jRhee & Katgert 


1987 


, 'West" 1989b' a'; 'Richstone et al.lll992l:|Plionis'l994') also 



reported both, alignment between neighbour clusters, as well 
as alignments with other groups embedded in the same su - 
percluster. For instance, using 48 superclusters lWesO (|l989al ) 
found a tendency for groups to be aligned on scales as large 
as 60 Mpc h"^ More recentlv lCodlowski fc Fliiil (|2010l ) anal- 



ysed the orientation of galaxy groups in the Local Super- 
cluster, finding strong correlations with the distribution of 
neighbouring groups up to scales of about 20 Mpc. Using 
the fourth Release of the Sloan Digita l Sky Survey (SDSS, 
I York fc et allboOOll . lWang et al.l (|2009l ') report several types 
of alignment signals between pairs of neighbouring galaxy 
groups and between groups and their surrounding galaxies. 

On the other hand, there is an extensive amount of work 
studying the presence of alignments between halo shapes 
and the large scale structure using numerical simulations 
( AraRon-Calvo et al. 2007; Hahn ct al. 2007; Bruni no et al. 
2007; Cuesta et al. 2008; P atiri ct a l. 2006; Basilakos et al. 
20061 ). Most of these studies perform statistics on the incli- 
nations of halo axes with respect to directions defined by the 
surrounding structure. This structure is characterised by the 
particle distribution through various topological signatures 
such as filaments, walls or voids. While these studies agree 
on the presence of alignments within the ACDM context, 
the magnitude of the alignment effect depends strongly on 
the definition of the surrounding structure. 

In this paper we report measurements of alignments 
between galaxy group shapes and the surrounding galaxy 
distribution in both, real galaxy groups obtained from the 
seven th data release of the Sloa n Digital Sky Survey (SDSS- 
DR7, lAbazaiian fc et al.ll2009l ) and from a numerical simu- 
lation. This is carrie d out following a methodology similar 
to those developed in lPaz et al.l |20o3), where the two-point 
correlation function was used to quantify the alignment be- 
tween galaxy angular momenta and the large scale struc- 
ture in both, numerical simulations and real data. In the 
present work we perform similar measurements, by using 
the "anisotropic" cross-correlation functions between galaxy 
groups and individual galaxies, in the directions parallel and 
perpendicular to the major axis of the group shape. Due 
to we are only interested in the alignment between galaxy 
groups and the surrounding large scale structure, we fo- 
cus our analysis on observed anisotropics over the two halo 
term of the two-point cross-correlation function. The advan- 
tage of this approach is the independence on any particular 
characterisation of the surrounding structure. Furthermore, 
this modified two-point cross-correlation function allows a 
robust comparison between numerical simulations and ob- 
ser vational data. A similar procedure t o the one presented 
by I Paz et al.l (120081 1 ■ was also applied bv lFaltenbacher et al.l 
l|2009l ). 



2 ACDM PREDICTIONS FOR HALO 
SHAPE-LSS ALIGNMENTS 

In this section we use a numerical simulation to study the 
anisotropy of the halo-particle cross-correlation function. 
The anisotropy signal around halos is characterised by tak- 
ing into account the relative orientation of a particle posi- 
tion with respect to the shape of the central halo. In order 
to mimic the observational effects on the anisotropy signal, 
we also analyse the distribution of galaxies around galaxy 
groups on a SDSS mock catalogue by means of the projected 
group-galaxy cross-correlation function. 

2.1 Numerical simulation 

Throughout this work we use a coUisionless numerical sim- 
ulation covering a periodic volume of SOO"' {h^^Mpcf . The 
initial conditions at redshift ~ 50 were calculated assum- 
ing a spatially flat low-density Universe, with a matter and 
vacuum density parameters Sim = 1 — = 0.258, Hub- 
ble constant Ho ~ 71.9 km s~^ Mpc~^, and normalisation 
parameter erg = 0.796. The resulting particle resolution is 
6.67 X 10^° h-'^ Mq. The run was performed using the 
second version of the GADGET code developed bv lSpringell 
(2005), with a gravitational soft ening of 0.03/i~^Mpc chosen 
by following [Power et al.l (|2003l ). 

The identification of particle clumps was carried out 
by means of a standard friends-of-friends algorithm with a 
percolation length given by Z = 0.17 where V is the 
mean number density of DM particles. In order to obtain 
reliable measurements of halo shapes (as is described below) , 
we only consider halos with at least 30 particle members. 
Consequently the resulting halo sample has a minimum mass 
threshold of ~ 2 x lO^'^ Mq. Nevertheless we are mainly 
interested in the observed galaxy group mass range, that is 
between 10^^ to 10^^ Mq (see figure [3}. Therefore our 
results can be associated with masses ranging from small 
galaxy groups to galax y clusters. 



Several authors ( 


Warren et al.l 19921: Thomas et al. 


19981: iHoDkins et al. 


20051: lAUgood et al. 


2006; 


Kasun & EvrardI 20051: iLau et al.1 I2OI0I: Paz et al.1 


2006, 



and references therein) have analysed the properties of 
halo shapes using the best-fitting ellipsoid to the spatial 
distribution of halo members. Following this standard 
method, we calculate the shape tensor for each dark-matter 
halo using the positions of each of its particle members. 
This can be written as a symmetric matrix, 

Uj^{i/Nh)"^Xc,acj, (1) 

where Xai is the i*** component of the displacement vec- 
tor of a particle a relative to the centre of mass, and Nh 
is the number of particles in the halo. The matrix eigen- 
values correspond to the square of the axis (a, b, c where 
a > 6 > c) of the characteristic ellipsoid that best describes 
the spatial distribution of the halo members. Our alignment 
analysis is based on the eigenvector directions (a, b and c). 
A fixed b/a = 1 with an arbitrary value of c/b, corresponds 
to perfect oblate ellipsoids. On the other hand, systems with 
fixed c/b = 1 are perfect prolate ellipsoids. A system with 
b/a < c/b is associated to a general triaxial ellipsoid with 
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Figure 1. Left panel: Spatial Halo— Dark-matter correlation function of low mass halos in the numerical simulation. The solid, long 
dashed and short dashed lines show the correlation functions between halos and neighbour particles along the directions of the major, 
intermediate, and minor axes of the halo shape ellipsoid, respectively. The dotted line shows the results for the isotropic correlation 
obtained by using all the tracers regardless of the direction (i.e. the classical halo-particle cross-correlation function). The dotted-dashed 
lines correspond to the correlation function between halos and neighbour particles in the direction perpendicular to the halo shape major 
axis. The dashed vertical line represents the transition scale between the 1— and 2— halo terms of the cross-correlation function. Shaded 
areas indicate the 1— cr standard deviation around each measurement. The lower left subpanel shows ratios between the directional and 
isotropic correlation functions shown in the top-left panel. Right panel: same as left panel, for a high mass DM halo subsample obtained 
from the simulation. 



prolate tendency, while the opposite case, b/a > c/b, corre- 
sponds to a predominantly oblate ellipsoid. 

Since we will apply a statistical method to observations, 
we also construct a synthetic SDSS catalogue in order to 
test the effects introduced by observational biases. By ap- 
plying the SDSS selection function on our sim ulation box 
we p roduce a mock catalogue (see for instance, I Cole et al] 
Il998t ). For simplicity, and because our goal is merely to 
match the geometry and depth of the catalogue, we com- 
pute the radial selection function of simulation particles by 
randomly assigning them an absolute magnitude. We use dif- 
ferent Schechter luminosity functions to apply SDSS r-band 
magnitudes to the particles. For particles outside groups we 
use the parameters « = 0.0149 M, = -20.44 a = -1.05 
l|Blanton et al.l [20031 ) . whereas for particles inside groups, 
we use different Sche chter parameters, whic h depend on the 
halo mass following (jZandivarez et al]|2006l ). The flux-limit 
selection of the mock catalogue is then obtained by impos- 
ing the same upper apparent magnitude threshold which 
affects the observational sample (in our case we consider 
galaxies with m < 17.5, in both real and mock catalogues). 
In order to determine the angular selection function of the 
survey, we divided the sphere in ~ 7, 700, 000 pixels using 
the SDSSPix software. This pixelisation scheme was specif- 
ically designed for the SDSS geometry (see for instance 
ISwanson et al.l |2008| ) . We then compute an initial binary 
mask by setting to 1 all pixels with at least one galaxy in- 
side, and the remaining pixels. Given that the pixel size 



is quite smaller than the mean inter-galaxy separation, we 
smooth the initial mask by averaging each pixel over its 
7x7 adjacent neighbours. We compound the angular selec- 
tion function as the set of all pixels with a value that exceeds 
a given threshold. Finally, we construct a SDSS-DR7 mock 
catalogue from the simulation box, by means of radial and 
angular selection functions. In addition, a mock group sam- 
ple is produced from the mock catalogue, by applying the 
same percolation algorithm used for real groups (see section 
I3.1|l . In our mock catalogue, the shapes of groups will be 
calculated in the same way as in the observational data, i.e. 
using only the projected positions of group members on the 
sky. 



2.2 Three-dimensional Results: anisotropic two 
point cross-correlation function 

We use the spatial halo-particle cross-correlation function 
^(r) to study the relation between the anisotropy presented 
by the distribution of matter and the orientation of halo 
shapes. This cross-correlation function measures the prob- 
ability excess to find a dark matter particle (tracer) in the 
volume clement dV at a distance r from a given halo (centre) 
(Peebles 1980). A standard method to measure ^(r) in a sim- 
ulation box consists in counting halo-particle pairs in a given 
distance bin, and normalising this quantity by th e expected 
numb er of pairs for a homogeneous distribution. |Paz et al] 
(|2008l ) perform these counts restricting the pairs to fixed 
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Figure 2. Global ratios between three dimensional correlation 
functions along the directions parallel and perpendicular to the 
central halo major axis {£,\\a/£.±a)- This ratio is computed by 
counting pairs with separations on the 2-halo scale regime. The 
horizontal dashed line corresponds to a no anisotropy behaviour. 



ranges in the angle between their relative position vector and 
a given direction. In our case, since we estimate the halo- 
particle cross correlation function, we restrict halo-particle 
pair counts to different ranges of the angle subtended by 
its relative position vector and the central halo shape axes 
(i.e. a, 6 or c). We will refer to the classical cross-correlation 
function, computed using all pairs regardless of their angle 
as the "isotropic" cross-correlation function (^iao). We define 
the "parallel" cross-correlation functions as those computed 
using all pairs subtending an angle with respect to a given 
halo shape axis smaller than a given threshold 6\ — 60° (we 
justify this value later in this section). These functions will 
be denoted and ^yc, for the corresponding axis a, b 

and c. 

We will also measure the dependence of the ampli- 
tude and shape of the parallel cross-correlation functions 
on halo mass (Mhaio). For this purpose, the halos in the 
simulation are split into five samples with different mass 
ranges. Figure [1] shows in dotted lines the isotropic halo- 
particle cross-correlation function ^iso for low mass systems 
(Mhaio ^ 2.5 X lO^^M0h~^, left panel) and high mass sys- 
tems (Mhaio > 1-3 X lO"M0h"\ right panel). The short 
dashed, long dashed and solid lines correspond to the cross- 
correlation functions estimated using particles in the direc- 
tion parallel to the minor, intermediate, and major axes of 
the halo shape, respectively, as described above. The bot- 
tom panels show the ratio between the directional cross- 
correlation functions and the isotropic function, with the 
same line styles as in the upper panels. Shaded areas corre- 
spond to Jacknife error bands {1-a uncertainties), estimated 
using the jacknife technique with a total of 100 subsamples 
for both, the numerical simulation and the SDSS data. Our 



jacknife errors estimations are stable between 20 to 150 sub- 
samples, without error underestimations. We will apply this 
error estimate throughout this paper. 

The halo-particle cross-correlation functions in Fig- 
ure [T] show the characteristic 1— and 2— halo terms re- 
ported in_several_2iwiou^^ observa- 
tional jZehavi fc et al.ll2004l : lYang et al1l2005h and numer - 
ical simulations (see for instance iHavashi fc Whi"t3 l2008l l. 
The dashed vertical line represents the transition scale be- 
tween 1— and 2— halo terms (where the correlation function 
has a local minimum in the first derivative). 

In Figure [1] it can be seen that the parallel cross- 
correlation functions along the minor and intermediate halo 
shape axes (^y^ and ^\\b), show a lower amplitude than the 
isotropic classical cross-correlation function (^iso), which in 
turn is smaller than the cross-correlation function measured 
along the halo major axis direction (Cya)- This measurement 
can be interpreted as an excess of surrounding large scale 
structure along the direction of halo major axis. The decre- 
ment of correlation power shown by ^||;, and ^yj,, could sig- 
nal a decrement of the particle density in the surrounding 
matter distribution in the plane perpendicular to the cen- 
tral halo major axis. The drop in the correlation amplitude 
with respect to the isotropic case observed for the ^yi, func- 
tion is smaller than that observed for the ^yc. This results 
reflects the triaxial nature of the structure surrounding a 
typical halo. In spite of this quantitative difference both, 
and exhibit a decrement rather than an increment, 
that could be understood as a weak preference for a triaxial 
over a prolate configuration for the surrounding structure. 
For instance, if would showed an increment with respect 
to the isotropic function in the same way as ^y^ does, we 
would have interpreted instead an oblate tendency for the 
surrounding structure. It is also important to remark that 
^l|b becomes consistent with the isotropic function ^iso at 
large enough scales. 

Given the prolate tendency shown by the cross- 
correlation function analysis, as well as the similar ampli- 
tudes of ^yi, and ^iso, from this point on, we will only study 
the anisotropy of the matter distribution around dark mat- 
ter halos using directions defined by the halo shape major 
axis. In addition to this, we compute the "perpendicular" 
cross-correlation function by adding pairs lying an angle 62 
or lower from the plane perpendicular to the halo major axis 
direction. This function is denoted by ^±a- We choose the 
threshold angles, for the parallel and perpendicular cases 
mentioned above {61 and 62, respectively), so that the vol- 
umes in each case are the same. This can be achieved by 
setting, sin(S2) = 1 — cos(Si) — x, and choosing a value for 
the threshold parameter x- When different values of x < 0.5 
are chosen, the angles 9i and 62 tend to acquire more acute 
values, while maintaining equal volumes. When adopting dif- 
ferent values of x we found that the ratios between parallel 
and perpendicular correlations, ^\\d/^±d (where d represents 
either axis, a, b or c) remain roughly constant, and the sig- 
nificance increases for larger volumes, up to x = 0.5. When 
the volumes for the parallel and perpendicular correlation 
are not the same, ^±d and ^y^j show different degrees of 
departure from the isotropic correlation; in particular, the 
correlation corresponding to the largest volume tends to be 
more similar to the isotropic case. Consequently, we have se- 
lected X ~ 0.5, which implies angles 81 — 60° and 62 = 30°. 
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This choice ensures that every halo-particle pair contributes 
to one of the two volumes. 

The results for the perpendicular cross-correlation func- 
tion are shown as dot-dashed lines in Figure [T] As can be 
seen, its amplitude is quite similar to the parallel cross- 
correlation function along the minor axis. In fact, over most 
of the 2-halo regime these two correlation functions are in- 
distinguishable according to the jacknife error bands. Addi- 
tionally, there is also an important difference between the 
parallel and perpendicular cross-correlation functions (with 
respect to the major axis of the halo shape, a); their rel- 
ative amplitudes show highly statistically significant differ- 
ences. The transition between the 1— and 2— halo terms can 
be clearly seen in the cross-correlation function quotients. 
Given the halo shape defin ition (Equation [T]) and the as- 
pherical geometry of halos l|Paz et al.l 120061 '). it is expected 
that on scales within the 1— halo term, the parallel and per- 
pendicular cross-correlation functions will present a signif- 
icant difference, since they will simply mirror the intrinsic 
halo shapes. A more relevant result is that coming from 
the comparison between the cross-correlation functions over 
the 2— halo term regime, which shows an approximately con- 
stant behaviour and only loses statistical significance at very 
large separations. Figure [T] also shows that the anisotropy 
of the cross-correlation function for the high mass sample 
(right panel) is higher than for the low mass sample. We 
have analysed five samples with halo masses ranging from 
IO^^Mq/i"^ to 2 X lO"Af0/i~\ confirming that the 2-halo 
term alignment signal increases with the halo mass. 

As is described above, the anisotropy signal behaviour is 
approximately constant over the 2— halo clustering regime. 
Therefore, it is possible to characterise a median alignment 
signal over this range by counting all pairs with separations 
in the 2— halo regime, in order to estimate a global quotient 
between correlation functions along the directions parallel 
and perpendicular to the halo centre major axis. These re- 
sults are shown on Figure [2] for five samples constructed us- 
ing different halo centre mass intervals. The median mass for 
each sample is shown on the abscissas, whereas the ordinate 
shows the global ratio (,\\a/S,±a- Error bars are estimated us- 
ing the jacknife resampling technique. As can be seen in the 
figure, there is an increment of the degree of anisotropy as 
the mass increases. This is expected in the hierarchical halo 
formation scenario since more massive halos have more re- 
cent formation epochs and therefore they are more likely to 
be aligned with the present large scale matter distribution. 

2.3 Projected anisotropic cross-correlation 
function and resolution effects 

The analysis of observational data from redshift surveys is 
mainly limited by three important effects, (i) The measure- 
ment of galaxy group shapes needs to be done in projection 
on the plane of the sky due to the presence of galaxy pe- 
culiar velocities which affect distance estimations based on 
redshift measurements. In order to asses the reliability of the 
estimates of projected group shape orientations, it is neces- 
sary to reproduce this effect in the simulation, and to study 
whether the alignment signals can still be detected in this 
case. The group shape tensor will therefore be calculated 
using only two directions along the plane of the sky. Due to 
the flux-limited nature of galaxy catalogues, (ii) only a few 



group members are available to estimate the shape tensor, 
and (iii) the number of tracers is diluted as the distance in- 
creases (both effects become stronger at higher redshifts). 
The decrement of galaxy group members (ii) tends to pro- 
duce deviations of the estimated group shape orientation 
with respect to the real direction. This is a shot noise type 
effect whic h also biases th e measured shapes towards oblate 
spheroids (Paz et al."200d) . The dilution of the number den- 
sity of galaxies with redshift can be taken into account with 
a proper modelling of the selection function of the catalogue 
via, for instance, constructing mock catalogues. We must 
estimate the influence of all these effects before proceeding 
further with the analysis of the real data. 

We start by projecting the simulated volume onto the 
Cartesian plane x — y in the numerical simulation. This plane 
resembles the "plane of the sky" of the observational data. 
The third Cartesian coordinate z plus the z-component of 
the peculiar velocity, is used to estimate a radial velocity 
along the "line of sight" . This last quantity mimics the red- 
shift measured in observations. The effects of low number 
statistics on the measured group shapes can be taken into 
account by selecting at random a small number of particle 
members in order to compute the projected shapes. Even 
though galaxies in groups do not necess a rily fo llow the dark 
matter halo density proflle, |Paz et al] (|2006l ) have shown 
that this procedure reproduces the biases present in the ob- 
servational estimates of group shapes due to projection and 
shot noise effects, that is, the shapes of simulated dark mat- 
ter halos measured this way exhibit indistinguishable differ- 
ences with shapes obtained from galaxy groups in both, real 
and mock catalogues. 

We now perform the equivalent of the analysis pre- 
sented in section [2?2] using anisotropic projected two-point 
cross-correlation functions. In an analogous manner to the 
method implemented on the three-dimensional case, we sep- 
arate halo-particle pairs according to the angle subtended 
with respect to the apparent halo major axis. Pairs with an- 
gles greater than and less than a threshold angle of 9p — 
45°, deflne the parallel (a;||) and perpendicular (ajj_) pro- 
jected cross-correlation functions respectively. Notice that 
this choice of Op is analogous to that of di and 62 for the 
three dimensional case, since this separates the plane of the 
sky in two equal areas. The quotients ijj^^/ijJj_, for smaller 
fractions of the area of the sky are roughly constant with a 
decreasing statistical signiflcance (due to a lower number of 
counts). This behaviour is analogous to the observed in the 
three dimensional signal. As a means to avoid a large amount 
of foreground and background structures we only consider 
neighbours in distance slices of 50 Mpc h~^ around the halo 
centres. Differences between the correlation amplitudes of 
and uj± will then be interpreted as statistical excesses of 
the surrounding large scale structure along the direction of 
the central group projected major axis. This will represent 
a lower limit of the real anisotropy around the group major 
axis given the averaging of the orientation of the 3D shape 
major axis with respect to the line of sight. 

With the purpose of studying the dependence of the 
anisotropy with the mass, we will measure the cross- 
correlation functions selecting halo centres separated into 
three mass ranges corresponding to the terciles of the Sloan 
group mass distribution (see section [3711 for the description 
of the group sample) . This choice allows a closer comparison 
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Figure 3. Distribution of group masses for different samples 
selected according to minimum thresholds on number of galaxy 
members. The lighter shaded histogram shows the results for the 
main galaxy group sample. Distributions of galaxy group subsam- 
ples with more than 15 and 20 galaxy members are shown with 
darker shaded areas as indicated in the figure key. 



with the observational data, and ensures the same number of 
galaxy groups in each mass range in the observational anal- 
ysis (see section r3.2|) . On the other hand, it is well known 
that the estimate of a halo shape depends on the number 
of members used to compute it l|Paz et al.l |2006| ). There- 
fore, it is expected that this will affect the alignment signal 
obtained from the ratios between anisotropic and isotropic 
cross-correlation functions. In order to analyse the extent of 
this effect we compute the halo shapes using only 10, 15 and 
20 random particles from each halo. Furthermore, to better 
mimic observations we impose that the distributions of halo 
meiss resemble the real galaxy group distributions with 10, 
15 and 20 members, respectively (see Figure 

Figure |4] shows the anisotropic projected halo-particle 
cross-correlation functions for the third tercile on the mass of 
groups with at least 10 members (Mhaio ^ l-7x IQ^^MqIT^) . 
The halo major axis direction was determined using 10 par- 
ticles selected at random. Solid, dotted and dashed lines 
correspond to parallel, isotropic and perpendicular cross- 
correlation functions to major axis, respectively. The lower 
panels show the quotients between anisotropic and isotropic 
cross-correlation functions. By comparing this figure with 
right panel of Figure [T] it can be seen that the alignment 
signal decreases by approximately a 75% due to shot noise 
and projection effects. It can also be seen that the quality 
of the signal decreases, lowering the statistical significance 
by some degree. However, the detection of the alignments 
remains strong. The error bars decrease significantly when 
20 particle members are used in the estimate of the major 
axis direction. 

So far, we have studied the combined effect of projection 
and shot noise; a third effect mentioned at the beginning of 
this section, the flux-limited nature of observational cata- 
logues, needs to be analysed. To this end, we use artificial 
data sets that resembles the real SDSS-DR7 catalogue (see 
Section 2.1) and the real galaxy group catalogue. 



R Scale [Upc h"'] 

Figure 4. Anisotropic projected correlation functions for mas- 
sive halos (third tercile) with a major axis direction determined 
using only 10 randomly selected particles. The mass distribution 
was forced to resemble the corresponding real galaxy group mass 
distributions with 10 members. 



Based on this mock data we compute again the parallel 
(a;|l) and perpendicular {(jJ±) projected group-galaxy cross- 
correlation functions, by applying the same procedure used 
on the observations, as will be described in the following 
sections (see Section I3.2|l . The results are shown in Figure 
[5] Line types are as in Figure |4] The first three panels (A, 
B and C) show the projected anisotropic cross-correlation 
functions for different mass terciles, for mock groups with at 
least 10 galaxy members. In the lower right panel of Figure[5] 
(panel D) we show the results for mock groups in the highest 
mass tercile with at least 20 group members. 

We now compare the results from projected and three 
dimensional correlation functions with the aim to assess the 
degree of loss of the alignment signal of section 12.21 The 
analysis of the expected biases allows us to conclude that the 
alignment signal obtained in section is still detectable on 
projected measurements of the two point cross-correlation 
function. When only considering the effect of projection we 
find, for scales around 10 Mpc h~^, differences between cor- 
relation quotients in the directions parallel and perpendicu- 
lar to the halo orientation decrease by a 50%, with respect 
to the corresponding three-dimensional ratio (Figure[T]right 
panel). When the mass distribution of observational galaxy 
groups is reproduced (Figure [Sjl , and also the discreteness 
effect on the determination of group shapes is taken into 
account (Figure |4} , the anisotropy signal diminishes even 
more, another 50% at scales around 10 Mpc h~^, that is a 
75% with respect to the quotients shown in Figure [T] Fi- 
nally, when the depth and geometry of the observational 
survey are taken into account, as well as the galaxy group 
identification is mimic, the error bars increase significantly. 
However, the anisotropy signal remains detectable with a 
high statistical significance over a wide range of distances. 
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Figure 5. Projected group— galaxy correlation functions for three group mass ranges extracted from a mock SDSS catalogue. The 
dashed lines correspond to the correlation function between mock groups and mock neighbour galaxies in the direction perpendicular to 
the group shape major axis, as seen projected onto the sky. The solid lines show the results when using tracers along the direction of the 
projected major axis, and the dotted lines show the corresponding results when all neighbours are taken into account. The first three 
panels. A, B and C, show the results for the lowest, intermediate and highest mass terciles of the mock galaxy group sample selected 
with at least 10 mock galaxy members per group. The lower right panel, labelled D, shows the results for the high mass tercile with at 
least 20 members (M ^ 2.6 X 10^^h~^MQ). The ratios between the correlation functions along different directions and the correlation 
function obtained using all the neighbours are shown in the lower subpanels. 



3 OBSERVED CORRELATION ALIGNMENTS 
ON SDSS REAL DATA 

In this section we analyse the presence of galaxy group 
shape alignments with the large scale structure as traced by 
real galaxies. To this end, measurements of anisotropic pro- 
jected cross-correlation functions were performed searching 
for an alignment signal in the 1— halo and 2— halo clustering 
regimes. At the end of this Section we provide a comparison 
between observational and numerical simulation results. 



3.1 The galaxy group sample 

We use the main galaxy sample from the 7th Data 
Release of the Sloan Digital Sky Survey (SDSS-DR7, 
'Abazaiia n fc et al.ll2009l ') which includes roughly ^ 900, 000 
galaxies with redshift measurements to z ~ 0.3 and an upper 
apparent magnitude limit of 17.77 in the r band. The galaxy 
groups used in this work have been identified in this sam- 
ple us ing the same procedure as in iMerchan fc Zandivarej 
(2OO5I). The method consists in using a friend-of-friend al- 
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gorithm similar to the one developed by iHuchra fc Gelled 
l| 19821 '). The algorithm links pairs of galaxies which sat- 
isfy Dij ^ D(iR(z) and Vij ^ VoR[z) where Dij is the pro- 
jected distance and Vij, the line-of-sight velocity difference. 
The scaling factor R{z) oc n^'' is introduced in order to take 
into account the decrement of the observed galaxy number 
density, Ug, with redshift due to the apparent magnitude 
limit cutoff (see eq. 5 in Ifuchra & GcUcr 1982). We have 
adopted a transverse linking length Do corresponding to an 
overdensity of 5p/p — 80, a line-of-sight linking length of 
Vb — 200 km s~^ and a fiducial distance of 10 Mpc h~^. Our 
main group sample comprises all galaxy groups with at least 
10 gala xy members. Group vir i al mas s es are estimated fol- 
lowing iMerchan fc Zandivarej l|2005l ): iLimber fc Mathewd 
(ll96C 



M : 



G ' 



where a is the three dimensi onal velocity dispersion, is 
the virial radius (see eq. 6 in lMerchan fc Zandivareal2005l ) 
and G is the gravitational constant. The velocity dispersion 
a is estimated using its observational counterpart, the line- 
of-sight velocity dispersion a^, a — \/3(t„. 

Figure[3]shows the mass distribution for the main group 
sample and for two additional subsamples with at least 15 
and 20 galaxy members. In order to analyse the dependence 
of any alignments with mass, we divide the ordered mass 
distribution in three parts, each containing a third of the 
sample population. The halo mass distribution shifts and 
stretches to higher masses when a higher minimum number 
of galaxy members is imposed. 



Observed alignments on the projected 
cross-correlation function 



3.2 



In section IzizI the correlation function measured in the sim- 
ulation was estimated by counting halo-particle pairs and 
normalising by the expected number of pairs for a homo- 
geneous distribution. Due to the complex selection function 
of the SDSS survey, it is necessary to produce a random 
distribution of points with the same selection function in 
order to normalise the group-galaxy pair counts. This nor- 
malisation can be done in several wa ys; in this work we use 
two estimators, the classic estimator (jPavis fc PeebledlTgsi ) 
uj = DD/DR — 1, where DD and DR are the numbers 
of group-galaxy and group-rando m pair counts re s pectiv ely, 
and a symmetric version of the iLandv fc Szalavl ()l993l l es- 
timator; in all cases we find negligible differences between 
these two estima tors. We perform all the analysis using the 
iLandv fc Szalavl estimator. This computation is performed 
as tJ = {DD-DR-RD+RR) /RR, where in addition to DD 
and DR, we need to calculate RR and RD, the numbers of 
random centre-random tracer and random centre-galaxy pair 
counts respectively. For each sample of groups we compute 
a random centre catalogue that mimic its redshift distribu- 
tion. The number of random centres is about 100 times the 
number of groups in each sample. The homogeneous sample 
of tracers was computed from the expected numerical den- 
sity of galaxies at a given redshift, with a magnitude below 
the limit of the survey, and a Schechter luminosity distri- 
bution with parameters (f)t = 0.0149 , M, = —20.44 ,a = 



-1.05 (jBlanton et all 120031 ). These random catalogues con- 
tain 2 X 10^ random points. The angular positions are re- 
stricted to the angular mask described in section [231 

In order to estimate the galaxy group anisotropic cross- 
correlation function we follow the same definitions given in 
section [2.31 (i.e. cj||, uJiso)- We therefore count DD and 
DR pairs depending on their relative orientation with the 
major projected axis of the group shape tensor. These counts 
are binned over a wide range of scales covering from 0.1 to 
50 Mpc h"'^. In analogy with the estimate of the projected 
cross-correlation function in simulations (section 12. 3|) . we 
only consider galaxy neighbours with relative comoving dis- 
tance < 50 Mpc h~^ around the group centre. 

In Figure[S]we show the resulting cross-correlation func- 
tion measurements. Line types are as in Figure O The first 
three panels (A, B and C) show the projected anisotropic 
cross-correlation functions for different mass terciles, for 
groups with at least 10 galaxy members. In order to study 
the effect of low number statistics, we repeat the same anal- 
ysis for samples with more than 15 and 20 galaxy members. 
For simplicity only the results from the highest mass tercile 
are shown in panel D of Figure (6] for groups with at least 
20 galaxy members. 

The first three panels in Figure |S1 show a similar be- 
haviour to that seen in the numerical simulation, where the 
anisotropy signal increases with mass, particularly for the 2- 
halo regime. In other words, significant statistical differences 
are seen between the and ujj_ projected cross-correlation 
functions, which become more prominent as the mass of the 
centre group increases. The number of galaxy groups used 
to compute the results showed in panel D is approximately a 
third of the total number of groups. Consequently, the cor- 
relation function errors increase. Nevertheless, the signal is 
still significant. A striking similarity of the results on this 
figure with the corresponding mock results can be observed 
by direct comparison with Figure [S] 

The results shown in this section confirm the presence of 
alignments in the large-scale structure as traced by galaxies 
and the direction of the major axis of galaxy groups, in qual- 
itative agreement with the expectations from ACDM model. 
This alignment signal persist over the two halo regime, up 
to scales of 20 Mpc , pointing an alignment between the 
group inner galaxy distribution and the surrounding large 
scale structure. 

It should also be mentioned that the covariance ma- 
trix of the correlation function has non-zero off-diagonal 
elements. In order to reduce the effects of covariance be- 
tween measurements at different scales, we estimated global 
ratios between cross-correlation functions in the directions 
perpendicular and parallel to the group major axes on two 
different wide scale ranges, corresponding to the 1— and 
2— halo terms, characterised by projected distances lower 
and greater than a transition scale. This procedure allows to 
diminish the effect of error correlations on the anisotropy, 
and also allows to study the alignments on these two dif- 
ferent scale regimes. The scale threshold depends on mass, 
varying from 0.5 to 3 Mpc h~^ over the mass range consid- 
ered here. 

Figure [7] shows the 1— and 2— halo term ratios as a 
function of the average subsample mass (left and right pan- 
els, respectively). In the figure the open symbols represent 
the results obtained from real observational data, whereas 
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Figure 6. SDSS projected group— galaxy correlation functions for three group mass ranges. The dashed lines correspond to the correlation 
function between groups and neighbour galaxies in the direction perpendicular to the group shape major axis, as seen projected onto 
the sky. The solid lines show the results when using tracers along the direction of the projected major axis, and the dotted lines show 
the corresponding results when all neighbours are taken into account. The first three panels, A, B and C, show the results for the lowest, 
intermediate and highest mass terciles of the galaxy group sample selected with at least 10 galaxy members per group. The lower right 
panel, labelled D, shows the results for the high mass tercile with at least 20 galaxy members (M ^ 2.6 X lO^^h^^MQ). The ratios 
between the correlation functions along different directions and the correlation function obtained using all the neighbours are shown in 
the lower subpanels. 



filled symbols correspond to mock data described in section 
12.31 As shown in the figure key, the different open (filled) 
symbols represent measurements obtained from real (mock) 
group samples with at least 10, 15 and 20 galaxy group mem- 
bers. This figure can be used to analyse the dependence of 
the anisotropy on the group mass for the 1— and 2— halo 
regimes in real and simulated data. 

The 1— halo term ratios obtained from the numerical 
simulation reflect a slight increase in the asphericity of halos 
as their mass increases, which is qualitatively consistent with 



previous results on the projected shapes of galaxy groups 
iPaz et al.ll|2006l ). The observations do not show a clear ten- 
dency, although error bars do not allow a clear distinction 
between mock and real data, indicating that the use of the 
shape tensor is more adequate for the purposes of studying 
the variation of the intrinsic shape of groups as a function of 
mass. In most cases the observations are in agreement with 
the numerical simulation results. 

On the 2— halo regime, the ratios shows a marginal in- 
crease with the group mass, both in the numerical Simula- 
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tion and observational results, which do not show a system- 
atic dependence on the minimum number of group members 
used. Within the uncertainties, the observational results are 
in agreement with mock data. It can be seen for all mass 
samples, regardless of the number of group members, a sig- 
nificant degree of anisotropy in the correlation ratios. For 
the more massive sample this anisotropy is distinguishable 
from isotropy at a level above 3 standard errors. 



4 CONCLUSIONS 

We have used a technique developed by iPaz et all (120081 1 
to characterise the alignment between the shape of gravi- 
tationally bound systems and the large scale structure on 
a ACDM numerical simulation and observed galaxy groups 
taken from the Sloan Digital Sky Survey, Data Release 7. We 
have detected a statistically significant alignment of galaxy 
groups with the surrounding large scale structure traced by 
SDSS spectroscopic galaxies. This result, based on an up- 
to-date sample of spectroscopic galaxy groups, supports the 
well known theoretical alignment of dark matter halo shapes 
with the circumventing large scale structure. 

The analysis on the ACDM numerical simulation shows 
that the direction of maximum elongation of dark matter 
halos is aligned with the surrounding matter distribution 
on scales of up to 30 Mpc h~^. The resulting halo-particle 
cross-correlation function shows the characteristic transition 
between the 1— and 2— halo terms regimes which is also re- 
flected in the anisotropy signal. Given that the 1— halo term 
corresponds to the inner matter distribution of the halos, 
the alignment on this range of scales was expected. More in- 
teresting to consider is the signal observed over the 2— halo 
term regime, which could be related with the halo forma- 
tion history. It has been established that halo formation oc- 
curs by anisotropic mergers with lower mass halos. These 
accretion events occur along directions traced to some ex- 
tent by the large scale structure. Therefore, the existence of 
a remnant correlation between halo shape orientations and 
the large scale matter distribution is expected. This con- 
stitutes a fossil record of the halo formation process. Due 
to the hierarchical nature of halo formation, an increasing 
anisotropy signal as the halo mass increases is also expected, 
as more massive halos finished assembling their mass more 
recently and therefore they are more likely to bear marks 
of more recent important formation processes. This mecha- 
nism also explains the sha pe-mass re l ation observed in the 
SDSS galaxy group sample [Paz et al.l l|2006l 'l. 



Our measurements of this alignment on the two dif- 
ferent halo model regimes are qualitatively consistent with 
previous results on dark matter halo alignmen t s in ACDM 



numerical studies ( 


Aragon-Calvo et al.' '2007': 'Hahn et al. 


2007: Brunino et al 




2007: Cucsta et al. 2008: Patiri et al. 


2009: iBasilakos et al 


EoodV The most extended method- 



ology along the literature, consists on performing statistics 
over different ranges in the relative angle defined by the halo 
major axis and the direction to larger structures d efined via 
partic ular geometrical criteria. For instance, Basilakos et al.l 
l|2006r ). found that cluster sized halos are aligned with their 
parent supercluster, in agreement with our results. The main 
asset of our approach is that it allows a robust comparison 
between numerical simulations and observational data, in- 



dependently of any particular characterisation of the sur- 
rounding structure. Given that the correlation function can 
be easily extended from three dimensions to projected data, 
this comparison arises naturally. 

Regarding the results in three dimensions, we have 
demonstrated that cross-correlation functions with tracers 
in the directions parallel to the minor and intermediate halo 
shape axes (^yc and £,\\b), show lower correlation amplitudes 
than the cross-correlation function measured along the ma- 
jor axis up to scales of 30 Mpc h~^. We also found 
important differences between the parallel and perpendic- 
ular cross-correlation functions (with respect to the major 
halo axis a). The detected correlation anisotropy signal was 
interpreted as a decrement of the typical particle density in 
the surrounding matter distribution in the plane perpendic- 
ular to the central halo major axis. We also showed that the 
alignment signal is higher as the central halo sample mass 
increases. 

In order to confront simulations and observations, we 
have analysed the effect of projections, group mass distri- 
bution and shot noise using the numerical simulation. The 
alignment signal remains detectable, with an important sta- 
tistical significance level over a wide range of distances. 

The key result of the present study is the observa- 
tional detection and quantification of a large scale anisotropy 
around galaxy groups. In agreement with full numerical 
simulations and mock catalogue results, our measurements 
of the projected two point cross-correlation function for 
SDSS groups and galaxies, exhibit an important degree of 
anisotropy. This signal is observed as a difference in the 
correlation amplitude estimated in the directions parallel 
and perpendicular to the group projected major axis. As 
expected, the amplitude of this feature increases with the 
average mass of the groups. The observed anisotropy is de- 
tected for group masses greater than ~ 6 x IO^'^Mq, and its 
statistical significance increases with the group mass. 

We average the alignment effect over ranges of scales 
corresponding to the 1— and 2— halo terms in the cross- 
correlation function. The 1— halo term results show similar 
anisotropics in the numerical simulation and the SDSS, and 
are consistent with t he known increa se of asphericity as the 
group mass increases (|Paz et al.ll2006l ) , a trend that is clearly 
detected in the numerical simulation and only mildly present 
in the observational sample of groups; this simply indicates 
that direct measurements of group shape (e.g. b/a) are more 
suited to detect this trend. 

The 2— halo regime shows an excellent agreement be- 
tween the mock model and observations. We find that the 
alignment of the shape of groups with the large-scale struc- 
ture increases with the halo mass reflecting the fact that 
higher mass groups have had, on average, more recent 
merger activity and therefore, younger dynamical ages. This 
result is independent of the number of group members used, 
indicating its lack of sensitivity to the shot-noise effect that 
biases group shape estimates (i.e. it affects the shape more 
strongly than the direction of the group major axis). The 
excellent quantitative agreement between observations and 
the ACDM model results provides new direct evidence of 
the adequacy of the gravitational instability to describe the 
large-scale structure formation of our Universe. 
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Figure 7. Global ratios between correlation functions along the directions parallel and perpendicular to the central group projected 
major axis (wy/wj,). Two different ratios are computed by counting pairs split into the 1- and 2-halo clustering regimes (left and right 
panels respectively) . The open symbols represent the results obtained from real observational data. Filled symbols show the corresponding 
results obtained from the SDSS mock catalogue. The horizontal dashed lines indicate the unit ratio (no anisotropy). 
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